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Abstract 



There is a striking convergence between Burgers turbulence and the con- 
tinuous spontaneous localization [CSL] model of quantum mechanics. In this 
paper, we exploit this analogy showing the similarities in the physics of these 
two apparently unrelated problems. It is hoped that the kind of analogy we 
introduce here may lead to important developments in both areas. 
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I. INTRODUCTION 



The study of nonlinear growth processes has recently become of central interest in physics, 
an example is provided by the KPZ equation |],|2|]. Originally devised to describe crystal 
growth, it has been since then applied to a wide range of systems, from bacterial growth ]3],|4j 
to directed polymers ||. One of the essential features of the KPZ equation is its nonlinear 
term introduced in order to account for lateral growth beyond the linear approximation such 
as described by the Edwards- Wilkinson model ||. 

The KPZ equation is given by 



(1) 



where /i(x, t) is the surface position, v is the viscosity and 0(x, t) is a space and time 
dependent noise. 

There is an interesting connection between the KPZ equation and the Burgers equation. 
The latter is a nonlinear diffusion equation for the velocity field of a fluid in N dimensions. 
The velocity is a gradient field, v = — Vh, and its equation is written as 
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where the external stirring force is given by f = — V0 and its correlation is 
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e is the energy injected into the fluid per unit time and unit mass and A is the length scale 
at which energy is injected. 

In terms of the random potential the correlation function is 



(0(x, t)0(x', £')) = eA 2 N5(t - t') exp 
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Through the Hopf-Cole transformation 0, v = — 2z/[VZ(x, i)]/Z(x, t), the Burgers equa- 
tion can be put into a linear form with a multiplicative noise 



dZ{x, t) 
dt 
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The above equation is a Schrodinger equation in the imaginary time. A similar equation 
for real time has been considered by several authors in a different context, such as the 
description of quantum open systems MM or the spontaneous collapse of the wave function 



[CSL model] in an attempt to solve the quantum measurement problem |]10| , p 



In a previous article ||12|| , we have presented an analogy between the CSL model and 
enhanced diffusion. We started from a comparison of our result for (x 2 ) in the beable 
interpretation of the CSL model with similar results obtained by Shlesinger et al |16| |. In 
the present paper, we show that this analogy extends much further by exploring a similarity 
between the CSL modified Schrodinger equation and the KPZ, Burgers equations. 
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The first section establishes an analogy between the KPZ and the CSL models, which 
provides a dictionary we use throughout the paper to move from one system to another 
and which allows us to introduce a Reynolds number for the CSL model. In section two 
we develop the stochastic picture underlying our beable interpretation of the CSL model 
and connect it with the Burgers equation through the velocity field. This enables us to 
introduce the intermittency corrections and a Fokker-Planck equation [FPE] corresponding 
to the stochastic equation for the CSL model in the random force approximation ||. An 
analysis of this system of equations allows us to introduce a time scale characterizing the 
dominance of enhanced diffusion over the standard (Wiener) diffusion. Some of the results 
are numerically illustrated with the sample paths of our stochastic differential equation. 
Section III presents our conclusions. 



II. A DICTIONARY FOR KPZ, BURGERS AND CSL 

The CSL model modifies the Schrodinger equation by introducing a multiplicative noise 
in the evolution of the wave function. In one dimension, the evolution equation in the 
Stratonovich form for a free particle of mass M is given by 

^&t^ = ~ X + VlJ dzw(z, t)G(x - z) | V>(x, t), (6) 

where w(z, t) is a white noise, so that (w(z, t)) = and (w(z, t)w(z', £')) = 5(z — z')5(t — t') 
and 



G(x -z) = ex P 



(x — z) 2 



—a- 



(7) 



characterizes the localization of the wave function. The diffusion constant v is equal to 
%/2M and 7 is related to the localization scale and the frequency of collapse A as 

7 = A \Ja-k /a. These parameters are chosen in such a way that the new evolution equation 
does not give different results from the usual Schrodinger unitary evolution for microscopic 
systems with few degrees of freedom, but when a macroscopic system is described there is 
a fast decay of the macroscopic linear superpositions which are quickly transformed into 
statistical mixtures fi~0| , [i~8| . 



The above modified Schrodinger equation is similar to eq.(g) in the imaginary time. From 
this analogy the noise term may be written as 4>(x,t) = lu^fy J dzw(z,t)G(x — z), which 
gives the following correlation function 



( ^2 

< <j>(x, t)(f>{x', t') >= Av 2 \5{t - t') exp[-a X ~ X) ]. (8) 



Comparing this equation with eq.(^), we identify the injection length scale A with y2/a 
and the injected energy e with 2u 2 a\ and can then easily obtain the Reynolds number for 
the CSL model by using the relation ||: Re = (eA 4 /z/ 3 ) 1//3 , which gives Re = (8A/ 'av) 1 ^ . 
For large Reynolds number we are in the domain of fully developed turbulence, which in the 
CSL case corresponds to a quantum system undergoing frequent collapses as A is large. 
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The analogy between turbulence and the CSL model also provides a clarification regard- 
ing the issue of energy non- conservation in the original model for the collapse of the wave 
function proposed by Ghirardi, Rimini and Weber [GRW] ||18|| . The celebrated Kolmogorov 
analysis of 1941 [K41] |17|| , identifies A as the scale at which energy e is injected into the 
system and a much smaller scale Id at which energy is dissipated. For distances r such that 
l(j < r < A, we have what is called the inertial regime, where energy is transferred to ever 
smaller length scales. The analogous result for the collapse model has that the energy is 
injected at the localization scale l/-\/a, which in turn is much larger than the atomic scale. 
In the framework of the collapse model, energy non- conservation coming from the collapse 
has always been discussed in terms of the magnitude of the effect itself and the constants 
were chosen in such a way as to make it a hardly observable effect at all. Here we stress, in 
analogy with turbulence, that the system does not gain energy indefinitely, for the energy 
injected in the collapse will be dissipated at the atomic scale. The disparity of scales is such 
that this fact was naturally unnoticed by the proponents of the collapse model. 

For completeness, let us end this section adding that it is in the inertial regime that 
Kolmogorov's scaling results 

(|v(x-r,t)-v(x,t)n~ri, (9) 

though it fails when p is larger than 3, a phenomenon known as intermittency, to which we 
shall return in the next section. 



III. DEVELOPING THE ANALOGY 

We have recently analyzed the CSL model from a microscopic point of view [|T3| . In order 



to do so we used Vink's treatment [pT9| , which shows that two alternative interpretations of 
quantum mechanics that treat position as a classical concept, the causal interpretation due 
to Bohm |2(| and the stochastic interpretation due to Nelson ||21]| , are actually particular 



cases of Bell's 'beable' interpretation f22|| . The beable interpretation is an attempt by Bell 
to treat physical quantities that exist independently of observation and therefore can be 
assigned well defined values. His approach used fermion number, a discrete quantity, and 
Vink extended it to any observable that takes discrete values on small scales. One starts 
from the equation for the probability density P m (t) on a given basis 

d t P m = Y, J mn, (10) 

n 

where the source matrix J mn is given by 

J mn (t) = 2Im{(^(t)\O m )(O m \H\O n }(O n m)}}. (11) 

From a stochastic point of view, the probability distribution of O m values satisfies the 
master equation 

dtP-m ^ / (T' mn P n T nm P m ^ (12) 

n 

where T mn dt is the transition probability for jumps from state n to state m. 
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To reconcile the quantum and stochastic views we equate fllOD and (|P^): 



Jmn (T mn P n T nm P m /j , 



(13) 



with T„,„ > and J„ 



" J 71. 



There is great freedom to find solutions of eq.(|K|). Bell chooses a particular one for 



T 



Jmn/ Pn Jmn ^* 
J mn ^ • 



Restricting the position of a particle in one dimension to the sites of a lattice, x = an, 
with n = 1,...,N and a the lattice distance, it follows from the discrete version of the 
Schrodinger equation that J mn is given by 



J r, 



Ma 



{[S(an)]'PJ n 



m—l 



[S(an)]'P n 5 n , m+ i}, 



(14) 



where use was made of the polar form of the wave function ip = Re^l^ and ip(x + a) is 
expanded up to first order in a. In the expression above [S'(an)]' = [S(an + a) — S(an)]/a. 
For forward movement, Bell's choice becomes 



T 



[S{an))' 
Ma 



n,m— 1 • 



(15) 



This term leads to transitions only between neighboring sites. However, we could also have 
added to T mn the solution of the homogeneous equation derived from eq.(|IB|): 



T° oc exp - 



m — n 



2aln{P m /P n 
4(m — n) 



/2a . 



(16) 



This extra term introduces transitions between more distant sites. 

Assuming a sufficiently small, we can approximate [ln(P m /P n )]/(m — n) by 
2a[R(an)]' / R(an) and arrive at the following Langevin equation for the particle position 
in the continuum limit 



x = v(x,t) + (f3aa ) 2 r](t), 



(17) 



where f3 is a free parameter, r](t) is a white noise, such that (f](t)) = 0, (7](t)i](t')) = 5(t — t') 
and 



v (x, t) 



((3aa 2 



dR(x,t) + 1 dS(x,t) 



R(x,t) dx M dx 



(18) 



Eq. (li~7|) coincides with Nelson's stochastic equation with flaa 2 = 2v. A similar equation 
was obtained in our analysis of the microscopic dynamics of the CSL model. This is so 
because the new terms in the source matrix J mn coming from the modified Schrodinger 
equation (|6|) do not contribute to the displacement dx |T^| . 
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The connection between the CSL (beable) equation and the Burgers equation is now 
straightforward. In imaginary time S(x,t) = and v = 2u(VR)/R, which, replacing R by 
Z, corresponds to the Hopf-Cole transformation. The velocity satisfies then the Burgers 
equation if we assume the relation between the random force and the random potential 
without the minus sign / = V0 (a point that was already noticed by Garbaczewski et al 

HI)- 

To simplify the following calculations, the Hamiltonian is set equal to zero as done in 



(19) 



nyi5|| . The solution of the modified Schrodinger equation is easily obtained 
ip(x,t) = exp(— At) exp y/j J dt' J dzG(x — z)w(z,t') 
By choosing the initial wave function as a single Gaussian 



V(x,0) 



(2nAx)* 



exp 



(x — (x)Y 



AAx 



(20) 



exp 



£ - x(x})\/AxAp - ti 2 /4 



Ax 



+ (p)x 



where Ax = (x 2 ) — (x) 2 and (x) is the mean value of x (the same is valid for p), the stochastic 
differential equation for position is 



x 



D s + 2v^ j jf* dt' J dzw{z, t') [-a(x - z)]G(x - z) 



2u7](t). 



(21) 



This equation describes the evolution of the position of a tracer in a turbulent medium. The 
first term on the right hand side describes a single free particle deterministic evolution and 
we choose as a short notation for the term derived from the initial wave function Ds = 
[2vVR s {x, 0)/R s (x, 0) + VS s (x, 0)/M] = [(x) - x] [v/Ax - yj AxAp - h 2 /A/MAx] + (p)/M. 
The two other terms describe the stochastic processes. The last term corresponds to a 
Brownian diffusion and the second one is responsible for the t 3 behavior of the mean square 
displacement, which is the same time dependence obtained by Richardson in his pioneering 
studies of turbulence [pi| . Notice that the coefficient of (x 2 ) corresponds to the injected 
energy 2a\v 2 as expected from hydro dynamical turbulence |23| and coincides with the 
result in the previous section. 

Moreover, the time dependence and nonlocal character of the second term on the r.h.s. 
of the equation above are in accordance with the concept of Levy walk as introduced by 



Shlesinger et al |T6| when studying the phenomenon of enhanced diffusion. The basic differ- 
ence between the more familiar Levy flight [[25] and Levy walk is that for the latter, although 
the walker visits all sites visited by the flight, the jumps do not occur instantaneously, but 
there may be a time delay before the next jump. By introducing time, Shlesinger et al ob- 
tained an integral transport equation involving a scaled memory which is nonlocal in space 
and time. Contrary to the infinite mean square displacement obtained in a Levy flight, the 
solution of such transport equation leads to a finite mean square displacement such as the 
one obtained by Richardson. 

From eq.fl2"ID, it follows that momentum satisfies 
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p = 2Mu^fy / dzw(z,t)[—a(x — z)]G(x — z). 



(22) 



The stochastic process for momentum derives from the introduction of the random potential 
in the modified Schrodinger equation, which in turn is responsible for the localization of the 
wave function. This process vanishes when the GRW parameters go to zero. 

The FPE corresponding to eq.fl2~l"D, which is a stochastic differential equation with a 

However 



colored noise introduced by v(t), is obtained as done by Zoller et al fl27lj28 
are in the inertial regime we use the random force approximation [ 2B|| : 
getting the following simplified equation for momentum 



The FPE is then 
dP(x,p, t) 
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(p) dP(x,p,t) 
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(24) 



A nice feature of our model is to obtain the above phase-space equation of evolution. 
Differently from Richardson |Tj| , we started from purely theoretical arguments and took into 
account the discontinuous nature of the particle velocity. Although eq. fl24"D is the same one as 
obtained by GRW [[jl|, we believe it is an oversimplification of the problem and the correct 
FPE to be used when studying the diffusion of a tracer should come from the complete 
equation for momentum, eq.(^). This is so because the random force approximation is 
deemed to be incorrect ||, since within it there is no intermittency. However, as we show 
below, there is a way out of this shortcoming which allows us to obtain the intermittency 
corrections even in the case of the random force approximation. 



A. Intermittency corrections 



Having eqs. fl21~D and ( p2|) we can now proceed to obtain the Mandelbrot intermittency 
corrections |29 to Richardson's law. In order to do so we replace a white noise in time by 
an affine one |30| called fractional Brownian noise 



(w(z, t) w{z', £')> = t A ~H{t - t') 5{z - z'), 



(25) 



which gives for the anomalous diffusion term [12 



(x 2 (t)) ~ t 



A- 1+3 



(26) 



This corresponds to one of the intermittency corrections obtained by Shlesinger et al [16 



provided we identify A — 1 with 3/i/(4 — /x), where \x = E — df, E being the Euclidean 
dimension and df, the fractal dimension. 

For the momentum variable the non-white noise gives 
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(p 2 ) - 1\ 



(27) 



which leads to the scaling relation obtained by Shlesinger et al [16| for the root-mean-square 
velocity. 

We now see that, within the random force approximation, the intermittency corrections 
can still be obtained as eq. (pleads to (r](t) r](t')) = t A ~ l 5(t - t'). 



B. Time scales 

,2\ 



Neglecting the deterministic term in eq. (PT|) , we find the mean square displacement (x ) 

(x 2 ) = 2ut + ~a\v 2 t 3 , (28) 
which has two contributions: the usual Brownian one, coming from the i](t) term, and the 



enhanced diffusion ||12|| , arising from the multiplicative noise in the modified Schrodinger 
equation. In fig.l we illustrate, with sample paths from eq. (|2T|) , the enhanced diffusion vis- 
a-vis the Brownian diffusion. We see that the numerical simulation exhibits the expected 
behavior of enhanced diffusion 

Comparison of the two terms in eq.(^) allows us to determine the time scale beyond 
which the enhanced diffusion dominates over the Brownian one 



tenh > J -J- (29) 

V aXu 



Using the GRW parameters a = 10 10 cm~ 2 , A(micro) = 10~ 16 s _1 , A(macro) 

= 10 7 s -1 , M(micro) = 10 -23 g, M(macro) = lg, we estimate the time scale, which is ap- 
proximately 2.4xl0 5 s, independently of the macro or microscopic nature of the system. 

A second time scale is given by the characteristic collapse time, which is of the order 
of A -1 . It does depend on the macro/microscopic nature of the system. Fig.2 displays the 
two time scales. The dashed line indicates the Schrodinger evolution, by which we mean 
the evolution of the wave function undisturbed by the collapse term. For the microscopic 
system enhanced diffusion manifests itself even before collapse occurs, opening an interesting 
window for experimental tests of this scenario, which would help to put stricter bounds on 
the parameters of the GRW/CSL models fS2[ . 



IV. CONCLUSIONS 

In this article we found further points of contact between turbulence and the CSL model, 
identifying the CSL counterparts of important quantities in turbulence, such as the Reynolds 
number and the injected energy. 

A microscopic stochastic picture for the CSL model allowed us to show that the velocity 
field of the tracer satisfies the Burgers equation. 

Finally a study of the different time scales governing enhanced and Brownian diffusion in- 
dicates that there is a region of parameter space of the CSL model amenable to experimental 
test, an investigation we plan to follow in a forthcoming publication. 
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FIG. 1. The curves were obtained using the algorithm for the numerical integration of stochastic 
equations with colored noise developed by Fox [31]. The horizontal axis corresponds to time in 
seconds. The full line represents the average value of x, which was obtained after 5000 realizations, 
and the dashed lines correspond to 4 sample paths. The parameters are those for a macroscopic 
particle and the time step is At = 10 4 s. 
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FIG. 2. The two lines give the time in seconds for the evolution of a microscopic system and a 
macroscopic one. The dashed lines correspond to the evolution of the wave function undisturbed 
by the collapse term. The thick line, which only appears for the microscopic system, corresponds 
to the time gap at which the CSL model may be tested. 



15 



